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Measurements of inclusive jet and dijet cross sections in photoproduction and deep inelastic 
scattering are presented. These measurements provide new tests of QCD, constrain the parton 
densities of the proton and the photon, and allow the search for new physics. Measurements 
of jet shapes are reported and used to test the differences between quark and gluon jets. 



1 Introduction 

The main sources of jets at HERA are deep inelastic scattering (DIS) and quasi-real photon- 
proton (7p) collisions (photoproduction). In DIS, a highly virtual photon (Q 2 3> 0, where Q 2 
is the virtuality of the exchanged photon) interacts with a parton from the proton. In hard 
photoproduction, a parton from the proton interacts with a quasi-real photon (Q 2 ~ 0) emitted 
by the electron beam. High-p^ jet production is sensitive to the parton densities in the proton, 
in particular the gluon density, and to the parton densities in the photon in "fp interactions. The 
measurements of jet cross sections provide constraints on the proton and photon parton densities, 
constitute new tests of QCD and allow the search for new physics. The internal structure of a jet 
is another tool to investigate QCD predictions and is sensitive to the differences between quark 
and gluon jets. During 1994 — 1997 HERA operated with positrons of energy E e = 27.5 GeV 
colliding with protons of energy E p = 820 GeV. 

2 Photoproduction at HERA 

In 7p interactions, the cross section for the production of jets at leading order (LO) in perturba- 
tive QCD (pQCD) is given by two contributions, namely those of resolved and direct processes. 
In resolved processes, a parton from the photon interacts with a parton from the proton, whereas 
in direct processes, the photon interacts as a point-like particle with a parton from the proton. 
The two contributions to the jet production cross section can be written as 

(^direct = J dn f 7 / e (y) fj/ p (x p ,n F ) da(jj -> jet jet) (1) 

and 



^resolved = J dn / 7 / e (y) fi/-f{x>y, fJ>jr) fj/ p {x p ,fip) da(ij -> jet jet), (2) 

where / 7 / e (y) is the flux of photons in the electron]] which is usually estimated by the Weizsacker- 
Williams approximation 111; fj/ p {x p ,n 2 F ) are the parton densities in the proton (determined from 
e.g. global fits!), x p is the fraction of the proton momentum taken by the parton and \ip is the 
factorisation scale; and da{^{i)j — > jet jet) is the subprocess cross section, calculable in pQCD. 
In the case of resolved processes, there is an additional ingredient: /i/ 7 (a? 7 , fJ-p) are the parton 
densities in the photon, for which there is only partial information. The integral is performed 
over the phase space represented by . 

^Talk given at the XXXIVth Rencontres de Moriond, QCD and High Energy Hadronic Interactions, Les Arcs, 
Savoie, France, March 20 th - 27 th , 1999. 

"The variable y is the fraction of the electron energy taken by the photon. 
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Figure 1: Measurement of the inclusive jet cross section da/dE^ in photoproduction. NLO QCD calculations 

are shown for comparison. 



3 Inclusive jet cross sections in photoproduction 

Inclusive jet cross sections in 7p interactions are sensitive to the underlying parton dynamics 
and to the parton densities in the photon. To test the available parametrisations of the photon 
parton densities, the experimental and theoretical uncertainties must be reduced as much as 
possible. Two approaches have been followed to reduce the experimental uncertainties due 
to the so-called "underlying event" !'§: e.g. to decrease the cone radius in iterative cone jet 
algorithms and/or to increase the jet transverse energy 13. In this way, a region of phase 
space was found where the next-to-leading order (NLO) QCD calculations are able to describe 
the dataEl. However, the uncertainties in the calculations for the iterative cone algorithm are 
still sizeable. The use of the kx cluster algorithm 13 reduces the theoretical uncertainties since 
it allows a transparent translation of the theoretical jet algorithm to the experimental set-up 
by avoiding the ambiguities related to the merging and overlapping of jets. In addition, a 
meaningful comparison between data and NNLO QCD calculations (when available) for the kx 
cluster algorithm will be possible. 

Inclusive jet cross sections have been measured i using the 1995 — 1997 ZEUS i data (which 
amounts to an integrated luminosity of C ~ 40 pb" 1 ) as a function of the jet transverse energy 
(E T et ). The jets have been searched for with the kx cluster algorithm. The measurements have 
been performed for jets of hadrons with E T et between 17 and 74 GeV and jet pseudorapidity 
{rf et ) between —0.75 and 2.5, and are given for the kinematic region defined by 0.2 < y < 0.85 
and Q 2 < 4 GeV 2 . 

Figure [l] shows the measured da/dE T et (black dots). In figures [l] and |2j, the systematic uncer- 
tainties not associated with the absolute energy scale of the jets have been added in quadrature 
to the statistical errors (thick error bars) and are shown as thin error bars. The shaded band 
represents the uncertainty on the energy scale of the jets. The data show a steep fall-off pf four 



usine 



orders of magnitude in the measured range. The curves are NLO QCD calculations L 
different parametrisatioiis of the photon structure function: GS96 £3 (solid line), GRV-HO 
(dashed line) and AFGli^ (dot-dashed line). The CTEQ4M&3 proton parton densities have been 
used in all cases. In the calculations shown here, the renormalisation and factorisation scales 
have been chosen equal to E T et and a s was calculated at 2-loops with A^L = 296 MeV. The 
NLO calculations give a reasonable description of the data. The right-hand side of figure |l] shows 
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Figure 2: Measurement of the high-mass dijet cross sections da/dM JJ and da/d\ cos#*| in photoproduction. NLO 

QCD calculations are shown for comparison. 



the fractional differences]] between the measured da/dE^ 1 and the NLO calculations. The NLO 
QCD calculations using the current knowledge of the photon structure are able to describe the 
data within the present experimental and theoretical uncertainties. 



4 High-mass dijet cross sections in photoproduction 

The dijet mass distribution M JJ is sensitive to the presence of new particles or resonances 
that decay into two jets. The distribution of the angle between the jet-jet axis and the beam 
direction in the dijet centre-of-mass system (cos 9*) reflects the underlying parton dynamics and 
is sensitive to the spin of the exchanged particle. New particles or resonances decaying into two 
jets may also be identified by deviations in the measured | cos 9* | distribution with respect to 
the QCD predictions. 

High-mass dijet cross sections have been measured using the 1995 — 1997 ZEUS data as a 
function of M JJ and | cos 9*\ in the kinematic region given by 0.2 < y < 0.85 and Q 2 < 4 GeV 2 . 
The measurements have been performed for M JJ > 47 GeV and | cos (9*1 < 0.8 using the kx 
cluster algorithm. 

Figure [2] shows the measured da/dM JJ (left-hand side) and da/d\ cos 9*\ (right-hand side). 
The data show a steep fall-off in M JJ of three orders of magnitude in the measured range. 
The measured da/d\ cos 9*\ rises as |cos#*| increases. The NLO QCD calculations give a 
reasonable description of the measured distributions. The predictions based on GRV-HO are 
closer in magnitude to the measured cross sections. No significant deviation between data and 
NLO calculations is observed in the measured range of M JJ and I cos a . 



5 Jet shapes in photoproduction 

The internal structure of the jets has been studied by means of the differential and the integrated 
jet shape. The differential jet shape (p(r)) is defined as the averaged fraction of the jet transverse 
energy which lies in the annulus of width Ar at a distance r from the jet axis, 



b The fractional differences are taken with respect to the NLO calculations based on GRV-HO. 
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Figure 3: Measured differential (left-hand side) and integrated (right-hand side) jet shapes in photoproduction. 
Calculations from the leading-logarithm parton-shower Monte Carlo PYTH1A are shown for comparison. 



p(r) - 1 T Mr-%,r + %) 

where Et(t — Ar/2,r + Ar/2) is the transverse energy within the given annulus and Nj e t s is 
the total number of jets in the sample. The integrated jet shape (ip(r)) is defined as 

QCD predicts that at sufficiently high transverse energy, the jet shape is driven by gluon emission 
off the primary parton and therefore calculable in pQCD. Gluon jets are predicted to be broader 
than quark jets due to the gluon-gluon coupling strength being larger than that of the quark- 
gluon. At LO, a jet consists of only one parton and thus has no structure. NLO QCD calculations 
give the lowest non-trivial order contribution to the jet structure. 

The differential and integrated jet shapes have been measured for jets defined with the cluster 
algorithm and required to have Ej? > 17 GeV and — 1 < if et < 2. The measurements have 
been performed using the 1995 ZEUS data (C ~ 6 pb _1 ) in the kinematic region given by 
0.2 < y < 0.85 and Q 2 < 1 GeV 2 . 

Figure [| shows the measured differential jet shape p(r) as a function of r in different regions of 
r] jet (left-hand side). The data (black dots) show that the jets become broader as if et increases. 
The solid histograms are the predictions of a leading-logarithm parton-shower Monte Carlo 
(MC) calculation using PYTHIAtH for resolved and direct processes. The calculations, which 
include initial- and final-state QCD radiation, give a good description of the data. The dashed 
(dot-dashed) histograms are the predictions of PYTHIA for samples of gluon (quark) jets. The 
measured jets are dominated by quarks for —1 < rf et < 0, and become increasingly more 
gluon-like as rf et increases. 

The quark and gluon content of the final-state jets has been studied in more detail by looking 
at the integrated jet shape for a fixed value of r (r = 0.3) as a function of rf et (right-hand 
side of figure ||). The measured jet shape decreases with r] jet , i.e. the jets become broader 
as rf et increases. Comparing the data to the model predictions for gluons (lower curve) and 
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Figure 4: Measured integrated jet shapes in DIS. On the right-hand side, calculations from the leading-logarithm 
parton-shower Monte Carlo LEPTO are shown for comparison. 



quarks (upper curve), it is observed that the data go from being dominated by quarks to being 
dominated by gluons as rf et increases. Thus, the broadening of the jets is consistent with an 
increasing fraction of gluon-initiated jets as rf et increases. 



6 Deep inelastic scattering at HERA 

Two processes contribute to the dijet cross section in neutral-current DIS at LO in pQCD, namely 
the QCD Compton (7*5 — > qg) and the Boson Gluon Fusion (BGF) (7*5 — > qq) processes. The 
cross section for jet production is given by 

ctdis = J dn fj/p(x p ,/j,%) da rj ^ et j et (x p ■ P, a s (/i 2 R ), /j? r , (5) 

where fj/ p (x p , n 2 F ) are the parton densities in the proton and dcr 7 * ;? _ > j e t j e t is the subprocess cross 
section. In the kinematic regimem studied, the BGF process represents the largest contribution 
and therefore the final state is expected to be dominated by quark- initiated jets. 



7 Jet shapes in DIS 

The integrated jet shape has been measured^ using the 1994 Hl@ data (£ ~ 2 pb" 1 ) for dijet 
events in the Breit frame. The jets have been searched for with the iterative cone and the kx 
cluster algorithms. The two jets with highest ^^Breit (transverse with respect to the direction 
of the virtual photon in the Breit frame) which satisfy E^^ Ieit > 5 GeV and — 1 < ??lab < ^ 
have been selected. The measurements are presented for different regions of 7/ B reit an< ^ ^TBreit 
in the kinematic region given by 10 < Q 2 < 120 GeV 2 and 2 • 10 -4 < xsj < 8 • 10 -3 , where xsj 
is the Bjorken x variable. 

Figure f| shows the measured integrated jet shape as a function of r in different regions of Threat 
and two regions of ^^Breit f° r the two jet algorithms. The data show that the jets get narrower 
as i^Breit increases and that the jets get broader as f/ B reit increases (i.e. towards the proton 
direction). The broadening of the jets is more evident at low Ej^ IfAt . At low E 3 T et BriAV the 
jets defined by the kx cluster algorithm are narrower than those defined by the iterative cone 
algorithm. At high EirBteit both algorithms produce very similar jets. 
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Figure 5: Measured dijet cross sections d 2 a/dQ 2 d£ and d 2 a/dQ 2 dxBj in DIS (left-hand side). NLO QCD calcu- 
lations are shown for comparison. On the right-hand side, the extracted gluon density in the proton is shown. 

The QCD based MC LEPTO fl predicts that gluon jets (dotted line in the right-hand side 
of figure H|) are broader than quark jets (dashed line) and that in the kinematic region of the 
measurements, the dijet sample is dominated by BGF and therefore by qq pairs in the final 
state. The comparison of the model predictions to the data shows that the measured jet shape 
is compatible with that of quark-initiated jets. 

8 Dijet cross sections in DIS 

Dijet cross sections in neutral-current DIS have been measured using the 1994 — 1997 HI 
data (C ~ 36 pb _1 ). The jets are defined by the far cluster algorithm in the Breit frame. 
The two jets with highest £^f* are required to have E^^^ > 5 GeV, £^ reit + E^^ mit > 17 

GeV and -1 < r/j^ B < 2.5. The measurements of d 2 a/dQ 2 d£, where £ = x B j(l + Mjj/Q 2 ), 
and d 2 a / dQ 2 dxBj are presented for the kinematic region given by 10 < Q 2 < 5000 GeV 2 and 
0.2 < y < 0.6. The data have been corrected for detector effects, and for initial- and final-state 
QED radiation effects. 

For the inclusive kj< cluster algorithm, the hadronisation effects are small and independent of 
Q 2 . Thus, the data have been compared to the NLO calculations performed using the program 

disentH. 

The left-hand side of figure || shows the measured dijet cross sections in DIS as a function of 
£ and XBj in different regions of Q 2 (black dots). The histograms are the NLO calculations 
from DISENT for two choices of the renormalisation scale: [x 2 R = Q 2 (solid histograms) and 
fi R = Q 2 + E^ (dashed histograms). The calculations based on different choices of fi R differ only 
at low Q 2 . The dependence of the calculations on the factorisation scale is small. The NLO 
predictions give a reasonable description of the data at all values of Q 2 for the choice /i R = Q 2 . 
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9 Extraction of the gluon density in the proton 



A direct determination of the gluon density has been performed via a NLO QCD fit to the 
measured dijet cross sections. Inclusive neutral-current DIS datac3 have been used to constrain 
the quark densities. The determination of the gluon density has been done using the dijet cross 
sections for Q 2 > 200 GeV 2 , where the hadronisation corrections are small and independent of 
Q 2 , £ and xsy The theoretical uncertainties are also small in this kinematic region. The value 
for the strong coupling constant has been taken from the world average, a s (M^) = 0.119±0.005. 
The parton densities have been extracted at a factorisation scale [i 2 F = 200 GeV 2 in the range 
0.01 < x < 0.1. 

The resulting quark densities are 5% higher at x = 0.01 than the results from global fits and 
they are in agreement at x = 0.1. The resulting gluon density together with its uncertainties 
is shown as a shaded band in the right-hand side of figure [|. This result is slightly higher than 
those from the global analyses, although compatible within the errors. It is in good agreement 
with the determination of the gluon density from the QCD analysis @ of F2 data at x = 0.01 
(dark solid line) by the HI Collaboration. 

10 Summary and conclusions 

Measurements of high-E'r inclusive jet and dijet cross sections in photoproduction have been 
presented. NLO QCD calculations give a reasonable description of the measured cross sections 
up to Ej? ~ 74 GeV. No significant deviation between data and calculations is observed up to 
dijet masses of M JJ ~ 140 GeV. 

The measurements of jet shapes in photoproduction and neutral-current DIS are consistent with 
the different parton contents of the final state expected in these two reactions: quark jets in 
DIS (7*5 — > qq is the dominant subprocess), and an increasing fraction of gluon jets as rf et 
increases in photoproduction (the resolved subprocess q^g v — * qg dominates over that of the 
direct 73 — > qq at large rf et ). 

A direct determination of the gluon density in the proton has been made via a QCD fit to the 
dijet and inclusive neutral-current DIS data. The resulting gluon density is compatible within 
errors with the results from global analyses, supporting the universality of the gluon density. 
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